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Abstract. New analytic solutions of the quadrupolc collective Bohr hamiltonian 
are proposed, exploiting an approximate separation of the fi and 7 variables to 
describe 7— soft prolate axial rotors. The model potential is a sum of two terms: 
a fi— dependent term taken either with a Coulomb-like or a Kratzcr-likc form, and 
a 7— dependent term taken as an harmonic oscillator. In particular it is possible 
to give a one parameter paradigm for a fi— soft, 7— soft axial rotor that can be 
applied, with a considerable agreement, to the spectrum of 234 U. 
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Figure 1. Contour plots of the Coulomb-like and Kratzer-like potentials in /3 
plus harmonic oscillator in 7. The coordinates are X = /3cos7 and Y = /3sin7. 
In the point fto = in the Coulomb-like case the potential tends to —00, while in 
the the Kratzer-like case tends to +00 (black area). The minimum in the right 
panel corresponds to (3q = 0.2 and D = 6. 



1. Introduction 

We have recently discussed PP a new analytic solution of the differential equation 
= Ety, where H is the Bohr hamiltonian with a 7— unstable class of potentials 
(called Coulomb- like and Kratzer-like potentials), in the very same spirit of the 
solution given in ref. [2] for the infinite square well. The hamiltonian discussed in 
the cited work displayed the dynamical symmetry SO(2,l)xSO(5) D SO(5)d SO(3) 
D SO(2). 

Now we would like to take a step further and exploit the approximate solution proposed 
in ref. [3] for the axially symmetric deformed rotor, replacing the infinite square well 
with our potentials for the (3 part of the problem. This is a mere exercise in the 
Coulomb-like case, because the peculiar behaviour of the (3 potential around zero is 
rather unrealistic in the nuclear case. The situation is more promising in the case 
of the Kratzer-like potential which displays a pocket with a deformed minimum and 
therefore, combined with an harmonic dependence on the 7 variable, may be taken as 
representative of (3— soft, 7— soft rotor. For the sake of illustration we give in fig. ^tw 
contour plots of the potential surfaces in {/3,7} plane resulting from the two cited (3 
potentials with an harmonic oscillator dependence on the 7 variable. From symmetry 
considerations 0] it is obvious that the 7 potential is restricted to be a function of 
cos 37. The choice of the harmonic oscillator is thus meaningful only around 7 = 0°, 
where it can be regarded as the first order Taylor approximation of the correct periodic 
even function. 

Following the procedure of ref. expanding the moments of inertia around 
7 = 0°, the Bohr hamiltonian may be recast in the following form 
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having expressed the cigcnsolutions as a product of a {[3, 7) part and of an angular 
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part (Wigner function) ^>((3, 7, 0j) = <p^-(/3, 7)2?m if(^j) an d having used the fact that 
within this approximation K is a good quantum number. As usual, reduced potential 
and energies, u = ^p-V and e = Wp-E, have been introduced, being B m the mass 
parameter of the Bohr hamiltonian. 
Under the assumptions that the potential is written as 

u(fi,i) = u(/3)+v(j) (2) 

and the 7 part has a minimum around 7 = 0°, the equation admits an approximate 
separation of variables with the factorized solution y^-(/3, 7) = Ci(/^) 7 7i<"(7) : 
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with e ~ 



e^ 7 ' and where (/3 2 ) is the average of (3 2 over 
In this note we will study the /3 part of the problem (in Section 2) with different 
choices of potentials, taking advantage (in Section 3) of Iachcllo's approximate solution 
of the equation in the 7— variable. Electromagnetic transition rates are evaluated (in 
Section 4). 

The peculiar feature of the present approach lies in the fact that the position of 
the ground state band and of all the (3— bands, as well as the corresponding moment 
of inertia along the bands, are fixed by only one parameter: an application of the 
solution with the Kratzer-like potential to the study of the spectrum of 234 U is given 
(in Section 5). 



2. (3— part of the problem 

With the substitution £l(/3) = xl(P)P~ 2 , equation J3J may be simplified to its 
standard form: 

r) 2 { (l I L ( L+1 > M 

W Xl{P) + |e<« - uC8) - f ^ \ \xl{P) = (5) 

The solution of this second order differential equation yields the eigenfunctions in 
the P variable and the eigenenergies e {/3 \ This is not the full solution of the problem, 
that requires also the solution of the equation in the 7— variable but represents the 
subset of the full spectrum corresponding to zero quanta in the 7 degree of freedom and 
lying lower in energy In this section we restrict our discussion to the (3— part, drawing 
a parallel between two cases corresponding to the two potentials, the Coulomb-like 
and the Kratzer-like potential, introduced in p^, leaving the considerations upon the 
7— part of the problem for the next section. 
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Figure 2. (3— spectrum of the Coulomb-like case (n 7 = 0). The vertical scale is 
expressed in units of the energy of the (v = 0, L = 2) state. Some B(E2) values 
are displayed in units of the lowest transition of the first band (numbers in italics). 
Notice the presence and energy of a threshold and the absence of degeneracy. 



2.1. Coulomb-like case 

Inserting the Coulomb-like potential u((3) = — A/(3 with A > and setting it in 
the variable x = with the further substitutions e = — e, k = and 

/i 2 = (| + £ 3 +1 ) ) w e obtain the Whittaker's equation 

d 2 , s \ 1 k 1/4- fi 2 ) , s 

^(*) + |- i + - + ^-|xL(x) = 0. (6) 

Its regular solution is the Whittaker function Mk )fl (x) that reads: 

X l(x) = Me- X ^x^ 1 F 1 (^+fi-k,l + 2 l i,x) (7) 

which is, in general, a multivalued function. The constant Af is determined from 
the normalization condition. We adopt usual conventions [3] and thus the function 
is analytic on the real axis. The hypergeometric series is an infinite series and to 
recover a good asymptotic behaviour we must require that it terminates, i.e. that it 
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Figure 3. Evolution of the ft— spectrum of the Kratzcr-likc rotor with the 
parameter B. The vertical scale is expressed in units of the energy of the 
(v = 0,L = 2) state. The B(E2) values, calculated with formula (1171 . between the 
lowest states are indicated in italics beside the corresponding downward arrows. 
Intraband transitions are not marked for the sake of simplicity. 



becomes a polynomial. This happens when the first argument is a negative integer, 
— i;, that must thence be regarded as an additive quantum number. This condition 
fixes unambiguously the (3 part of the spectrum: 

e 0*) = A*/4 

Qualitatively this spectrum (fig. [2J looks rather similar to the one studied in the 
case of 7— instability. The states are grouped into bands corresponding to different 
values of v; within each v all the angular momenta consistent with the constraints 
given by the 7— part (see Section 3) are present. Note however that, fixing to unity 
the energy difference of the first two states of the lowest band (eo,2 — £0,0)7 the energy 
of first 4 + state is 1.718, while the energy of the (v = l,L = 0) state (beginning of 
the second band) is 1.423. These values should be compared with the energy (1.35) 
of the corresponding degenerate states of the 7— unstable case of ref. PQ. Other pe- 
culiarities are the presence of the threshold, corresponding to an infinite L quantum 
number, that is located at 2.5616 and the absence of the degeneracy pattern found in 
the 7— unstable case. Only accidental degeneracies may occur, but they are quite rare 
(the lowest is between the L = 14, v = and L = 0, v = 7 states). 
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2.2. Kratzer-like case 

We consider now the case of the Kratzer-like potential. Inserting in eq. (J3J) the 
potential 

u(P) = -A/0 + B/13 2 (9) 

and setting x = e = — e, k = and fi 2 = (B + | + L ^ L § \ we obtain 

again the Whittaker's equation whose solutions may be written as in eq. JjJ and we 
may repeat the whole procedure of the previous section, the only major difference 
being the definition of /i. The spectrum assumes in this case the form 

A 2 /A 



t v,L 




(10) 



We recall that the two parameters used here, A and B, have an immediate translation 
into the position, /3o, and depth of the minimum of the potential, T>, being valid the 
following relations: A = 2/3qT> and B = (3 2 V . Consequently (3q = 2B/A and 
V = A 2 /(AB). 

In fig. |21 we show the excitation spectrum of the Kratzer-like potential in units of 
(eo.2 — £o.o)- This spectrum is clearly independent of A. We thus study its evolution 
with the parameter B. As expected, when B is very small we recover the Coulomb- 
like spectrum of fig. |3 while when B becomes very large we recover a typical 
rotational spectrum, as can be showed either numerically or analytically (for large 
B, (eo,L — £o,o) oc L(L + \)/{B 2 )). Among the two limits the ratio (relevant 
for a comparison with experimental data) varies between the two extremes, 1.718 and 
10/3. 



3. 7 part of the problem 

The approximate solution of the 7— part of the problem follows closely the derivation 
in 3 . In this reference the differential equation (0} in the 7 variable with an harmonic 
oscillator dependence upon 7 may be approximativcly solved by expanding all the sine 
functions in powers of 7. Thus the differential equation to be solved is 
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for eigenvalues, and 



(13) 
(14) 

for eigenfunctions in terms of Lagucrre polynomials, with n = (h 7 — | K |)/2 and 
n 7 = 0, 1, 2, 3, ... . A simple rule to express in a compact way the set of K quantum 
numbers encountered for a given n 7 is 



Vn-,,K(l) = 7 
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Figure 4. Kratzcr potential in j3 for 7 = in energy units of (to, 2 — £0,o)- The 
parameters (specified in the two equivalent sets on the figure) have been chosen in 
such a way that the lowest states of the first two bands (labeled by their quantum 
numbers v, L) are deeply bound in the well. 



The values of L are determined by K: if K = then L — 0,2, 4, 



while if K ^ 



then L = K, K + 1, K + 2, 



Levels in fig. [21 and [21 correspond to n 7 = 0, that 



implies K = 0, and therefore only even values of L are present. States labeled by the 
same set of quantum numbers, disregarding the sign of K, are degenerate in energy. 

4. Electromagnetic transitions 



The quadrupolc electromagnetic operator reads 

1 



T 



t(3 



Z>jg(fi)coB 7 - 



V2 
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(fi) ) sin 7 
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where t is a scale factor. Limiting ourselves to transitions within the n 7 = bands 
and exploiting the expansion of the trigonometric functions around 7 = 0, transition 
rates may be evaluated separating the integral in three parts: the angular integral is 
done analytically, exploiting the properties of the Wigner T> functions, the integral in 
7 reduces to an orthonormality condition, while the integral in (3 has been calculated 
numerically employing the explicit forms of the wavefunctions, 



B(E2,ViLi -> v f Lf) r 
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Figure 5. Reduced energies normalized to the energy difference between the first 
two states of the ground state band as a function of the parameter B. The first 
three states of the g.s. band (v = 0) as well as the first two of the /3— band (v = 1) 
are displayed. One can notice that the beginning of the second band lies at higher 
energies when the deformation of the system is increased and eventually tends to 
infinity in the rotor limit. 
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The calculated B(E2) values (normalized to the value for the transition from the first 
2+ state to the + ground state of the lowest band) are reported in fig. 0and|31 In 
the case of Kratzer potential, the B(E2) values tend to the Coulomb-like case values 
when the parameter B is small and tend to the rigid rotor values when B is large (for 
example in the case of the 4^" — > 2|J" transition we get the limiting value 10/7). 



5. Discussion 



The (/?, 7) potential with the Kratzer form in the (3 variable and the harmonic 
behaviour in 7 may be assumed as a representative of (3— soft, 7— soft rotor. For 
illustration we display in fig. 0] a section of a Kratzer potential along the line with 
7 = 0°. The parameters have been chosen to show how the well can be shaped in order 
to have a rotational behaviour in presence of a f3— softness. In fact the position of some 
of the lowest states of the first two bands is quite deep in the well: this means that the 
variation of the mean square deformation will not change much among these states (so 
that the approximate separation of variables in eq. JU is a posteriori justified). Note 
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Figure 6. Comparison between the calculated and experimental spectrum of 
234 U (positive parity states). The theoretical predictions are obtained according 
to CSJ and Git with thc choice B = 393.7415 and c = -41.5629 to define the 
(/3, 7) potential in JTJJ and 1121 . The ground state, (3 and 7 bands are displayed. 



that the present solution requires only one parameter to fix the relative position of all 
the states in the lowest band and the position of the so-called (3— bands (the parameter 
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A in the potential does not play any role in determining the scaled eigenvalues). The 
position of the 7— band is fixed by the additional parameter c in cq. I|13|) . As usual 
the physical values of the spectrum are then obtained by an overall scaling factor. 
As we move to large deformations and deeper minima by increasing the parameter 
B, the position of the j3— bands with respect to the ground state band moves higher 
in energy and a crossing occurs between the relative positions of states in different 
bands, as one can see in fig. 

We applied our form of the potential to describe the lowest rotational bands in 
234 U. A comparison between the lowest positive parity experimental energy levels of 
234 U [H] and the predictions of the model presented here are given in fig. The 
agreement is good in spite of the fact that all the relative positions of the states in 
the ground state band and (3— band have been fixed with only one parameter, B, that 
has been obtained fitting the energy of the 4+ level. 

Note that, at variance with the experiment that seems to display almost equal 
moments of inertia for the ground and f3— band, the theory gives a slightly larger 
moment for the (3— band. This can be understood from the expansion of the general 
formula l|10|) in powers of L(L + 1). In leading order one obtains, for the energy 
differences within a generic band (characterized by the quantum number v), the 
following expression 

A 2 /24 

e v ,L-e vfi = L(L+T)—= — - ' , (18) 

y/B + 9/4( + + 1/2 + v) s 

which leads to the mild dependence on v, discussed above. 

In summary we have discussed a new solution of the Bohr hamiltonian with the 
aim of describing a /3—soft, 7— soft axial rotor. We have given spectra and transition 
rates that may in principle encompass a broad range of situations and can be used to 
survey experimental data. We have discussed, as an example, the spectrum of 234 U, 
that may be considered as a "quasi-rigid" rotor and that is accurately described within 
the present formalism, making use of the solution with a Kratzcr-likc potential. 
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